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Abstract —High detection complexity is the main impediment 
in future Gigabit-wireless systems. However, a quantum-based 
detector is capable of simultaneously detecting hundreds of 
user signals by virtue of its inherent parallel nature. This in 
turn requires near-capacity quantum error correction codes 
for protecting the constituent qubits of the quantum detector 
against the undesirable environmental decoherence. In this quest, 
we appropriately adapt the conventional non-binary Extrinsic 
Information Transfer (EXIT) charts for quantum turbo codes by 
exploiting the intrinsic quantum-to-classical isomorphism. The 
EXIT chart analysis not only allows us to dispense with the 
time-consuming Monte-Carlo simulations, but also facilitates the 
design of near-capacity codes without resorting to the analysis 
of their distance spectra. We have demonstrated that our EXIT 
chart predictions are in line with the Monte-Carlo simulations 
results. We have also optimized the entanglement-assisted QTC 
using EXIT charts, which outperforms the existing distance 
spectra based QTCs. More explicitly, the performance of our 
optimized QTC is as close as 0.3 dB to the corresponding hashing 
bound. 

Index Terms —Quantum Error Correction, Turbo Codes, EXIT 
Charts, Near-Capacity Design. 


I. Introduction 

Multi-User Multiple-Input Multiple-Output (MU- 
MIMO) m, El and massive MIMO 0 schemes are 
promising candidates for the future generation Gigabit- 
wireless system. However, the corresponding detection 
complexity increases exponentially with the number of users 
and antennas, when aiming for approaching the optimum 
Maximum-Likelihood (ML) performance. An attractive 
solution to this exponentially escalating complexity problem 
is to perform the ML detection in the quantum domain, since 
quantum computing allows parallel evaluations of a function 
at a complexity cost that is equivalent to a single classical 
evaluation a, 0 . However, a quantum detector requires 
powerful Quantum Error Correction codes (QECC’s) for 
stabilizing and protecting the fragile constituent quantum bits 
(qubits) against the undesirable quantum decoherence, when 
they interact with the environment a, a. Furthermore, 
quantum-based wireless transmission is capable of supporting 
secure data dissemination a. a. where any ‘measurement’ 
or ‘observation’ by an eavesdropper will destroy the quantum 
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entanglement, hence intimating the parties concerned 11. 
However, this requires powerful QECC’s for the reliable 
transmission of qubits across the wireless communication 
channels. Hence, near-capacity QECC’s are the vital enabling 
technique for future generations of wireless systems, which 
are both reliable and secure, yet operate at an affordable 
detection complexity. 

Classical turbo codes operate almost arbitrarily close to the 
Shannon limit, which inspired researchers to achieve a com¬ 
parable near-capacity performance for quantum systems 0- 
11121 . In this quest, Poulin et al. developed the theory of 
Quantum Turbo Codes (QTCs) in 0, 0, based on the 
interleaved serial concatenation of Quantum Convolutional 
Codes (QCCs) |[T3l-lfT6i, and investigated their behaviour on 
a quantum depolarizing channeQ It was found in 0, IfTTl 
that the constituent QCCs cannot be simultaneously recursive 
and non-catastrophic. Since recursive nature of the inner code 
is essential for ensuring an unbounded minimum distance, 
while the non-catastrophic nature is required for achieving 
decoding convergence, the QTCs of 0, 0 had a bounded 
minimum distance. More explicitly, the design of Poulin et 
al. 0, |9j was based on non-recursive and non-catastrophic 
convolutional codes. Later, Wilde and Hseih iflOll extended the 
concept of pre-shared entanglement to QTCs, which facilitated 
the design of QTCs having an unbounded minimum distance. 
Wilde et al. also introduced the notion of extrinsic information 
to the iterative decoding of QTCs and investigated various 
code structures in i m. 

The search for the optimal components of a QTC has been 
so far confined to the analysis of the constituent QCC dis¬ 
tance spectra, followed by intensive Monte-Carlo simulations 
for determining the convergence threshold of the resultant 
QTC, as detailed in [0, QTj. While the distance spectrum 
dominates a turbo code’s performance in the Bit Error Rate 
(BER) floor region, it has a relatively insignificant impact 
on the convergence properties in the turbo-cliff region IT8l . 
Therefore, having a good distance spectrum does not guarantee 
having a near-capacity performance - in fact, often there is 
a trade-off between them. To circumvent this problem and 
to dispense with time-consuming Monte-Carlo simulations, in 
this contribution we extend the application of EXIT charts to 
the design of quantum turbo codes. 

1 A quantum channel can be used for modeling imperfections in quantum 
hardwares, namely faults resulting from quantum decoherence and quantum 
gates. Furthermore, a quantum channel can also model quantum-state flips 
imposed by the transmission medium, including free-space wireless channels 
and optical fiber links, when qubits are transmitted across these media. 
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More explicitly, 

• We have appropriately adapted the conventional non¬ 
binary EXIT chart based design approach to the family of 
quantum turbo codes based on the underlying quantum- 
to-classical isomorphism. Similar to the classical codes, 
our EXIT chart predictions are in line with the Monte- 
Carlo simulation results. 

• We have analyzed the behaviour of both an un-assisted 
(non-recursive) and of an entanglement-assisted (recur¬ 
sive) inner convolutional code using EXIT charts for 
demonstrating that, similar to their classical counter¬ 
parts, recursive inner quantum codes constitute families 
of QTCs having an unbounded minimum distance. 

• For the sake of approaching the achievable capacity, 
we have optimized the constituent inner and outer com¬ 
ponents of QTC using EXIT charts. In contrast to the 
distance spectra based QTCs of am whose performance 
is within 0.9 dB of the hashing bound, our optimized QTC 
operates within 0.3 dB of the capacity limit. However, 
our intention was not to carry out an exhaustive code 
search over the potentially excessive parameter-space, 
but instead to demonstrate how our EXIT-chart based 
approach may be involved for quantum codes. This new 
design-approach is expected to stimulate further interest 
in the EXIT-chart based near-capacity design of various 
concatenated quantum codes. 

This paper is organized as follows. Section QI] provides 
a rudimentary introduction to quantum stabilizer codes and 
QTCs. We will then present our proposed EXIT-chart based 
approach conceived for quantum turbo codes in Section Hill 
Our results will be discussed in Section HVl while our conclu¬ 
sions are offered in Section [V] 


II. Preliminaries 

The constituent convolutional codes of a QTC belong to the 
class of stabilizer codes ED, which are analogous to classical 
linear block codes. We will here briefly review the basics of 
stabilizer codes in order to highlight this relationship for the 
benefit of readers with background in classical channel coding. 
This will be followed by a brief discussion of QTCs. 


A. Stabilizer Codes 

Qubits collapse to classical bits upon measurement 0 , 
0 . This prevents us from directly applying classical er¬ 
ror correction techniques for reliable quantum transmission. 
Quantum error correction codes circumvent this problem by 
observing the error syndromes without reading the actual 
quantum information. Hence, quantum stabilizer codes invoke 
the syndrome decoding approach of classical linear block 
codes for estimating the errors incurred during transmission. 

Let us first recall some basic definitions 0 . 


Pauli Operators: The I, X, Y and Z Pauli operators are 
defined by the following matrices: 



where the X, Y and Z operators anti-commute with each 
other. 

Pauli Group: A single qubit Pauli group Q\ consists of all 
the Pauli matrices of Eq. 0} together with the multiplicative 
factors ±1 and ±i, i.e. we have: 


Qx = {±1, ±*I, ±X, ±*X, ±Y, ±*Y, ±Z, ±i.Z}. (2) 


The general Pauli group Q n is an n-fold tensor product of Q\. 

Depolarizing Channel: The depolarizing channel charac¬ 
terized by the probability p inflicts an n-tuple error V £ Q n 
on n qubits, where the i th qubit may experience either a bit 
flip (X), a phase flip (Z) or both (Y) with a probability of 
p/3. 

An [n, k] Quantum Stabilizer Code (QSC), constructed over 
a code space C, is defined by a set of (n — k) independent 
commuting n-tuple Pauli operators gi, for 1 < i < (n — k). 
The corresponding encoder then maps the information word 
(logical qubits) | ip) £ C 2 onto the codeword (physical qubits) 
\'f) £ C 2 , where C d denotes the c?-dimensional Hilbert 
space. More specifically, the corresponding stabilizer group % 
contains both g, and all the products of < 7 , for 1 < i < (n — k) 
and forms an abelian subgroup of Q n . A unique feature of 
these operators is that they do not change the state of valid 
codewords, while yielding an eigenvalue of —1 for corrupted 
states. Consequently, the eigenvalue is —1 if the n-tuple Pauli 
error V anti-commutes with the stabilizer //, and it is +1 if V 
commutes with //,. More explicitly, we have: 


9i 


1^) 



gfP = Vgi 
gfP = -Vgi, 


where V is an n-tuple Pauli error, | ip) £ C and | ip) =V\f>) is 
the received codeword. The resultant ±1 eigenvalue gives the 
corresponding error syndrome, which is 0 for an eigenvalue of 
+1 and 1 for an eigenvalue of —1. It must be mentioned here 
that Pauli errors which differ only by the stabilizer group have 
the same impact on all the codewords and therefore can be 
corrected by the same recovery operations. This gives quantum 
codes the intrinsic property of degeneracy 11201 . 

As detailed in ETTl . l22l . QSCs may be characterized in 
terms of an equivalent binary parity check matrix notation 
satisfying the commutativity constraint of stabilizers. This can 
be exploited for designing quantum codes with the aid of 
known classical codes. The (n — k) stabilizers of an [n, k] 
stabilizer code can be represented as a concatenation of a pair 
of (n — k) x n binary matrices H z and H^, resulting in the 
binary parity check matrix H as given below: 


H=[H 2 |H,.]. (4) 

More explicitly, each row of H corresponds to a stabilizer of 
TL, so that the i th column of H 2 and H x corresponds to the 
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i th qubit and a binary 1 at these locations represents a Z and 
X Pauli operator, respectively, in the corresponding stabilizer. 
Moreover, the commutativity requirement of stabilizers is 
transformed into the orthogonality of rows with respect to the 
symplectic product defined in [ 22 ], as follows: 

H Z H^ + H.Hj = 0. (5) 

Conversely, two classical linear codes H- and H, can be used 
to construct a quantum stabilizer code H of Eq. © if H- and 
H x meet the symplectic criterion of Eq. ©. 

In line with this discussion, a Pauli error operator P can 
be represented by the effective error P, which is a binary 
vector of length 2 n. More specifically, P is a concatenation 
of n bits for Z errors, followed by another n bits for X 
errors and the resultant syndrome is given by the symplectic 
product of H and P, which is equivalent to H [P x : P z } 7 . 
In other words, the Pauli-X operator is used for correcting 
Z errors, while the Pauli-Z operator is used for correcting X 
errors | 6 ]. Thus, the quantum-domain syndrome is equivalent 
to the classical-domain binary syndrome and a basic quantum- 
domain decoding procedure is similar to syndrome based 
decoding of the equivalent classical code Ii22l . However, due 
to the degenerate nature of quantum codes, quantum decoding 
aims for finding the most likely error coset, while the classical 
syndrome decoding finds the most likely error. 

B. Quantum Turbo Codes 

Analogous to classical Serially Concatenated (SC) turbo 
codes, QTCs are obtained from the interleaved serial concate¬ 
nation of QCCs, which belong to the class of stabilizer codes. 
However, it is more convenient to exploit the circuit-based 
representation of the constituent codes, rather than the con¬ 
ventional parity-check matrix based syndrome decoding 11231 . 
Before proceeding with the decoding algorithm, we will briefly 
review the circuit-based representation. This discussion is 
based on 0 . 

Let us consider an (n, k) classical linear block code con¬ 
structed over the code space C, which maps the information 
word ceFj onto the corresponding codeword c £ FJf. In the 
circuit-based representation, the code space C is defined as 
follows: 

C = {c={c: 0„_ fe )V}, (6) 

where V is an (n x n)-element invertible encoding matrix 
over F 2 . Similarly, for an [n, k] quantum stabilizer code, the 
quantum code space C is defined as: 

C = { $> = V(|V0 ® |0„- fc »}, (7) 

where V is an n-qubit Clifford transformatior@ and b/>)G<C 2 \ 
The corresponding binary encoding matrix V is a unique (2 n x 
2 n)-element matrix such that for any V & Q n we have 0: 

[VPV f ] = [P]V, (8) 

2 Clifford transformation V is a unitary transformation, which maps an 72- 
qubit Pauli group Q n onto itself under conjugation ED, i.e. 

= g n . 



Encoder Inverse Encoder 


Fig. 1. Quantum transmission model. 



Fig. 2. Circuit representation of the inverse encoder PV 1 = (L : S). 


where [P] = P and [.] denotes the effective Pauli group G n 
such that [P] differs from V by a multiplicative constant, i.e. 
we have [P] = P/{±1,±/}. The rows of V, denoted as V) 
for 1 < i < 2 n, are given by V, = [VZ,V'\ = [Zi]V for 
1 < i < n and V t = [VA'jVl] = [Xj]V for n < i < 2 n. Here 
Xi and Z, represents the Pauli X and Z operators acting on 
the i th qubit. Furthermore, any codeword in C is invariant by 
VZiVK for k < i < n, which therefore corresponds to the 
stabilizer generators gi of Eq. (©. More explicitly, the rows 
Vi, for k < i < n, constitute the (n — k) x 2 n parity check 
matrix H of Eq. ©, which meets the symplectic criterion of 
Eq. (©. 

At the decoder, the received codeword | ip) = V\i[)) is passed 
through the inverse encoder IP, which yields the corrupted 
transmitted information word C\ijj) and the associated syn¬ 
drome <S|0 n _/ r ), formulated as: 

v f p|$) = v f pv(|v>) ® |o„_ fc )) 

= (cm g) (S|0„_fc)), (9) 

where C denotes the error imposed on the logical qubits, while 
S represents the error inflicted on the remaining ( n—k ) qubits. 
This transmission process is summarized in Fig. [I] 

Since stabilizer codes are analogous to linear block codes, 
syndrome decoding is employed at the receiver to find the most 
likely error coset C given the syndrome S. This is efficiently 
achieved by exploiting the equivalent binary encoding matrix 
V of Eq. ©. which decomposes the effective n-qubit error 
imposed on the physical qubits P = [P] into the effective 
fc-qubit error inflicted on the logical qubits L = \C\ and the 
corresponding effective (n — fc)-qubit syndrome S = [5], as 
portrayed in Fig. [2] and mathematically represented below: 

PV- 1 = ( L : S). (10) 

More explicitly, P £ G n , L £ Gk and S £ G n -k- 
For an [n, k, m] QCC, the encoding matrix V is constructed 
from repeated uses of the seed transformation U shifted by 
n qubits, as shown in Fig. 6 of 0 . More specifically, U is 
the binary equivalent of an (n + m)-qubit symplectic matrix. 
Furthermore, Eq. ( ITOt may be modified as follows 0 : 

(P t : MJU- 1 = (M t _! : L t : S t ), (11) 

where t and (t — 1 ) denotes the current and previous time 
instants, respectively, while M is the effective m-qubit error 














4 


on the memory states. Furthermore, 2(n — fc)-element binary 
vector S of Eq. m and (ITTll can be decomposed into two 
components, yielding S = S x + S z , wher S x and S z are 
the X and Z components of the syndrome S, respectively. 
The (n — A:)-hinary error syndrome computed using the parity 
check matrix H only reveals S x but not S z 0. Therefore, 
those physical errors which only differ in S z do not have to 
be differentiated, since they correspond to the same logical 
error L and can be corrected by the same operations. These 
are the degenerate errors, which only differ by the stabilizer 
group as discussed in Section Hl-AI Consequently, a quantum 
turbo decoding algorithm aims for finding the most likely error 
coset acting on the logical qubits, i.e. L , which satisfies the 
syndrome S x . 

Similar to the classical turbo codes. Quantum turbo decod¬ 
ing invokes an iterative decoding algorithm at the receiver 
for exchanging extrinsic information ED. E3 between the 
pair of SC Soft-In Soft-Out (SISO) decoders, as shown in 
Fig. |3] These SISO decoders employ the degenerate decoding 
approach of 0. Let Pi and /., denote the error imposed 
on the physical and logical qubits, while Sf represents the 
syndrome sequence for the i th decoder. Furthermore, P“(.), 
P|(.) and P°(.) denote the a-priori, extrinsic and a-posteriori 
probabilities ll25l related to the i th decoder. Based on this 
notation, the turbo decoding process can be summarized as 
follows: 

• The inner SISO decoder uses the channel information 
Prh(I\ ), the a-priori information gleaned from the outer 
decoder Pf(Li) (initialized to be equiprobable for the 
first iteration) and the syndrome Sf to compute the 
extrinsic information P\{L\). 

• V\(Li) is passed through a quantum interleave^ (7r) to 
yield a-priori information for the outer decoder P'f(Pi)- 

• Based on the a-priori information P'UJ'Y) and on the 
syndrome Sf, the outer SISO decoder computes both 
the a-posteriori information Pf(J-'i) and the extrinsic 
information PKP 2 ). 

• PKP 2 ) is de-interleaved to obtain P“(Li), which is fed 
back to the inner SISO decoder. This iterative procedure 
continues until convergence is achieved or the maximum 
affordable number of iterations is reached. 

• Finally, a qubit-based MAP decision is made to determine 
the most likely error coset L 2 . 

III. Application of EXIT Charts to Quantum Turbo 

Codes 

In this section, we will extend the application of EXIT 
charts to the quantum domain, by appropriately adapting 
the conventional non-binary EXIT chart generation technique 
to the circuit-based quantum syndrome decoding approach. 
Some of the information presented in this section might seem 
redundant to the experts of classical channel coding theory. 
However, since EXIT charts are not widely known in the 

3 An .V-qubiL quantum interleaver is an .V-qubiL symplectic transformation, 
which randomly permutes the N qubits and also applies single-qubit sym¬ 
plectic transformations to the individual qubits (£)■ 



Sf Sf 


Fig. 3. Schematic of the quantum turbo decoder. P“(.), P((.) and P°(.) 
are the a-priori, extrinsic and a-posteriori probabilities related to the i th 
decoder; Pi and Li denote the error on the physical and logical qubits, 
while Sf represents the syndrome sequence for the 7 , decoder. 


quantum community, this introduction was necessary to make 
this treatise accessible to quantum researchers. 

A. EXIT Charts 

EXIT charts ff 8 l . |25l , li26l are capable of visualizing the 
convergence behaviour of iterative decoding schemes by ex¬ 
ploiting the input/output relations of the constituent decoders 
in terms of their average Mutual Information (MI) charac¬ 
teristics. They have been extensively employed for designing 
near-capacity classical codes |27)-|29)- Let us recall that the 
EXIT chart of a serially concatenated scheme visualizes the 
exchange of the following four MI terms: 

• average a-priori MI of the inner decoder, l\, 

• average a-priori MI of the outer decoder, l\, 

• average extrinsic MI of the inner decoder, and 

• average extrinsic MI of the outer decoder, Jj,. 

More specifically, l\ and /),, constitute the EXIT curve of 
the inner decoder, while l\ and l\ yield the EXIT curve 
of the outer decoder. The MI transfer characteristics of both 
the decoders are plotted in the same graph, with the x and y 
axes of the outer decoder swapped. The resultant EXIT chart 
quantifies the improvement in the mutual information as the 
iterations proceed, which can be viewed as a stair-case-shaped 
decoding trajectory. Having an open tunnel between the two 
EXIT curves ensures that the decoding trajectory reaches the 
( 1 , 3 /) point of perfect convergence. 

B. Quantum-to-Classical Isomorphism 

Before proceeding with the application of EXIT charts for 
quantum codes, let us elaborate on the quantum-to-classical 
isomorphism encapsulated in Eq. @, which forms the basis of 
our EXIT chart aided approach. As discussed in Section III-AI 
a Pauli error operator V experienced by an A'-qubit frame 
transmitted over a depolarizing channel can be modeled by an 
effective error-vector P, which is a binary vector of length 2 TV. 
The first TV bits of P denote Z errors, while the remaining TV 
bits represent X errors, as depicted in Fig. 0] More explicitly, 
an X error imposed on the I s * qubit will yield a 0 and a 1 
at the 1 st and (TV + 1 ) th index of P, respectively. Similarly, 
a Z error imposed on the 1 st qubit will give a 1 and a 0 
at the 1 st and (TV + \) th index of P, respectively, while 
a Y error on the 1 st qubit will result in a 1 at both the 
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Fig. 4. Effective error (P) corresponding to the error imposed on an N -qubit 
frame (P). 
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TABLE I 

Correlation between X and Z errors on th i th qubit in terms 

OF THE CORRESPONDING PROBABILITY OF OCCURRENCE. 


1 st as well as (N + l) th index of P. Since a depolarizing 
channel characterized by the probability p incurs X, Y and 
Z errors with an equal probability of p/3, the effective error- 
vector P reduces to two Binary Symmetric Channels (BSCs) 
having a crossover probability of 2p/3, where we have one 
channel for the Z errors and the other for the X errors. 
Hence, a quantum depolarizing channel has been considered 
analogous to a Binary Symmetric Channel (BSC) E2 . (30), 
whose capacity is given by: 

Cbsc = 1 — H 2 (2p/2>), (12) 

where H 2 is the binary entropy function. Using Eq. ©I, we can 
readily infer that the code rate Rq of an [n, k] QSC is related 
to the equivalent classical code rate Rc as follows ED, ED: 

Rc = 2 (l + Rq) ■ (13) 

Consequently, the corresponding quantum capacity is as fol¬ 
lows (22), ED: 

Cg c = l-H 2 (2p/3). (14) 

However, the two BSCs constituting a quantum depolarizing 
channel are not entirely independent. There is an inherent 
correlation between the X and Z errors 11221 . which is char¬ 
acterized in Table Q] This correlation is taken into account 
by the turbo decoder of Fig. [3] Alternatively, a quantum 
depolarization channel can also be considered equivalent to a 
4-ary symmetric channel. More explicitly, the i th and ( N+i) th 
index of P constitute the 4-ary symbol. The corresponding 
classical capacity is equivalent to the maximum rate achievable 
over each half of the 4-ary symmetric channel, as follows (22, 

ED: 

C 4 . ary = i[2 - H 2 (p) — plog 2 (3)]. (15) 

Therefore, using Eq. ( IT3l > the corresponding quantum capacity 
can be readily shown to be ED, ED: 

cfary = 1 ~ H 2 (p) - p log 2 (3), (16) 

which is known as the hashing bouncfl. 

4 Hashing bound determines the code rate at which a random quantum code 
facilitates reliable transmission for a particular depolarizing probability p tm 



Fig. 5. System model for generating the EXIT chart of the inner decoder. 


Recall that a quantum code is equivalent to a classical 
code through Eq. ©. More specifically, as mentioned in Sec- 
tion lll-Al the decoding of a quantum code is essentially carried 
out with the aid of the equivalent classical code by exploiting 
the additional property of degeneracy. Quantum codes employ 
syndrome decoding H23H , which yields information about the 
error-sequence rather than the information-sequence or coded 
qubits, hence avoiding the observation of the latter sequences, 
which would collapse them back to the classical domain. 

Since a depolarizing channel is analogous to the BSC and 
a QTC has an equivalent classical representation, we employ 
the EXIT chart technique to design near-capacity QTCs. The 
major difference between the EXIT charts conceived for the 
classical and quantum domains is that while the former models 
the a-priori information concerning the input bits of the inner 
encoder (and similarly the output bits of the outer encoder), 
the latter models the a-priori information concerning the 
corresponding error-sequence, i.e. the error-sequence related 
to the input qubits of the inner encoder L\ (and similarly error- 
sequence related to the output qubits of the outer encoder P 2 ). 
This will be dealt with further in the next section. 


C. EXIT Charts for Quantum Turbo Codes 

Similar to the classical EXIT charts, in our design we 
assume that the interleaver length is sufficiently high to ensure 
that 01, E3: 

• the a-priori values are fairly uncorrelated; and 

« the a-priori information has a Gaussian distribution. 

Fig. 0 shows the system model used for generating the EXIT 
chart of the inner decoder. Here, a quantum depolarizing 
channel having a depolarizing probability of p generates the 
error sequence Pi, which is passed through the inverse inner 
encoder Vj . This yields both the error imposed on the logical 
qubits L\ and the syndrome Sf according to Eq. ( I IQl i. The 
a-priori channel block then models the a-priori information 
P“(Li) such that the average mutual information between 
the actual error L\ and the a-priori probabilities Pj (L \) is 
given by I a (Ja) fl8) . (25) , (26) . More explicitly, we have 
Ia(L i) = I[L i, P“(Li)], where I denotes the average mutual 
information function. As discussed in Section 1III-BL the i th 
and (N + i) th bits of the effective error vector L\ can be 
visualized as 4-ary symbols. Consequently, similar to classical 
non-binary EXIT charts (32 , (33), the a-priori information is 
modeled using an independent Gaussian distribution with a 
mean of zero and variance of a\, assuming that the X and Z 
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errors constituting the 4-ary symbols are independent. Using 
the channel information P c h(Pi), syndrome Sf and the a- 
priori information, the inner SISO decoder yields the extrinsic 
information Pf(Li) based on the classic forward-backward 
recursive coefficients at and Pt as follows 0: 

• For a coded sequence of duration N, let 

Pi = [Pi,i, Pi, 2 , ■ ■ ■, Pi,t, ■ ■ ■, Pi,n] and L l = 

[Li,i,L lt 2 ,...,L 1>t ,...,L liN ], where P M G G n and 
Li, t G G k . More explicitly, P M = [Pl t , P [ 2 V ..., P{\ t ] 
and Li jt = [L\ t , L\ t ,..., L\ t ]. For the ease of 

clarification, we will ignore the first subscript, which 
represents the decoder, in the algorithm given below, i.e. 
we have P\ = P and L\ - L. Similarly, Sf = S x . 

• Let U = (Up : Um ), so that Up is the binary matrix 

formed by the first 2 n columns of U of Eq. (ITTIl . while 

Um is the binary matrix formed by the last 2 m columns 
of U. Therefore, we have: 

Pt = (Mt~i : L t : S t ) U P . (17) 

M t = (M t -1 : L t : S t ) U M ■ (18) 


Finally, the extrinsic average mutual information between L\ 
and Pf(Li) is computed, i.e. Ie(L\) = I[L\, Pf (L\)\. Since 
the equivalent classical capacity of a quantum channel is 
given by the capacity achievable over each half of the 4- 
ary symmetric channel as depicted in Eq. 03. Ie(Li) is the 
normalized mutual information of the 4-ary symbols, which 
can be computed based on ll33l . Il34l as: 


X;p;(Li (m) )i 0 g 2 pf(Li (m) ) ), 
m=0 J } 

(24) 

where E is the expectation (or time average) operator and 
L : ^ m> is the 4-ary m th hypothetical error imposed on the logical 
qubits. More explicitly, since error on each qubit is represented 
by an equivalent pair of classical bits, is a 4-ary classical 
symbol with m G {0, 3}. The process is repeated for a range of 
Ia(Li) G [0,1] values for obtaining the extrinsic information 
transfer characteristics at the depolarizing probability p. The 
resultant inner EXIT function T) of the specific inner decoder 
may be defined as follows: 


Ie(Li) = - 2 


• Let at (M t ) be the forward recursive coefficient, which 
is defined as follows: 


I E (L 1 )=T 1 (I A (L 1 ),p\, (25) 


a t (M t )4p(M t |S| t ), 

oc Y, P ° = A ) p ch (Pt) a t -1 (p ), (19) 

/X,A,£7 

where p G G m , A G G k and cr G G n - k , while 
a = <j x + a z , having a x = Sf. Furthermore, we have 
P t = (p : A : cr)Up and M t = (p : X : a) Um- The 
channel information P, :/ , (P t ) is computed assuming that 
each qubit is independently transmitted over a quantum 
depolarizing channel having a depolarizing probability of 
p, whose channel transition probabilities are given by 0: 

p.fpU-J i-P’ ifP f =I ( 20) 

ch ' * > \ p/3, if VI G {X, Z, Y}. ( 

• Let fit (M t ) be the backward recursive coefficient, which 
is defined as follows: 

Pt (M t ) — P (M t \Sf t ), 

oc Y pa ( L ‘ = A ) P ch (Pt+ 1 ) Pt+ 1 (M t +i ), (21) 

A ,<j 

where Pt+i = (M t : A : cr)Up and M t + 1 = 
(M t : X : a) U M - 

• Finally, we have the a-posteriori probability P °(L t ), 
which is given by: 

P°(L t )±P(L t \S x ), 

oc Y pa (L t )Pch(Pt)a t -i (p) Pt (M t ), (22) 

fi,cr 

where P t = (p : L t : a)Up and M t = (p : L t : a) Um- 

• Marginalized probabilities P °(L J t ) for j G {0, k — 1} are 
then computed from P")!-) ) and the a-priori information 
is removed in order to yield the extrinsic probabilities 
CD, i.e we have: 

ln[P e (^')] = ln[P°(^')] - ln[P a (L\)\. (23) 


which is dependent on the depolarizing probability p of the 
quantum channel. 

The system model used for generating the EXIT chart of 
the outer decoder is depicted in Fig. [6] As inferred from the 
figure, the EXIT curve of the outer decoder is independent of 
the channel’s output information. The a-priori information is 
generated by the a-priori channel based on P -2 (error on the 
physical qubits of the second decoder) and I A (P 2 ), which is 
the average Mutual Information (MI) between P 2 and P') () • 
Furthermore, as for the inner decoder, P 2 is passed through 
the inverse outer encoder U 2 -1 to compute Sf, which is fed 
to the outer SISO decoder to yield the extrinsic information 
P e 2 (P 2 ). Based on Eq. ( 1 1 9b and ( 12 1 1 . this may be formulated 
as follows 0: 

P°(P t )±P(P t \S x ), 

(X Y p ( p t) p ( L t = A H- 1 (p) Pt (Mt) , (26) 

where P t = (p : A : o)Up and M t = (p : X : a) Um- 
The resultant probabilities are marginalized and the a-priori 
information is removed similar to Eq. (l23l >. The average MI 
between P 2 and P^iPi) is then calculated using Eq. ( 1241 ). 
The resultant EXIT chart is characterized by the following MI 
transfer function: 

Ie(P 2 ) = T 2 [I a (P 2 )}, (27) 

where T 2 is the outer EXIT function, which is dependent 
on the specific outer decoder, but is independent of the 
depolarizing probability p. 

Finally, the MI transfer characteristics of both decoders 
characterized by Eq. d25l > and Eq. ( 127b are plotted in the same 
graph, with the x and y axes of the outer decoder swapped. 
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Fig. 6. System model for generating the EXIT chart of the outer decoder. 


PTOIR (Outer) - PTOIR (Inner) 



Fig. 7. The EXIT curves of a QTC parametrized by the increasing 
depolarizing probability p (rate-1/9 QTC having PTOIR as both the inner 
and outer components was used). 


IV. Results and Discussions 
A. Accuracy of EXIT Chart Predictions 

In order to verify the accuracy of our EXIT-chart based 
approach, we have analyzed the convergence behaviour of 
a rate-1/9 QTC, consisting of two identical rate-1/3 QCCs. 
More specifically, for both the inner and outer decoders, we 
have used the configuration termed as “PTOIR” in tm no, 
which is a non-catastrophic but quasi-recursive code. 

Our first aim was to predict the convergence threshold using 
EXIT charts, which would otherwise require time-consuming 
Word Error Rate/Qubit Error Rate (WER/QBER) simulations. 
Convergence threshold can be determined by finding the max¬ 
imum depolarizing probability p, which yields a marginally 
open EXIT tunnel between the EXIT curves of the inner and 
outer decoder; hence, facilitating an infinitesimally low QBER. 
Fig- 0 shows the EXIT curves for the inner and outer decoders, 
where the area under the EXIT curve of the inner decoder 
decreases upon increasing p. Eventually, the inner and outer 
curves crossover, when p is increased to p = 0.13. More 
explicitly, increasing p beyond 0.125, closes the EXIT tunnel. 
Hence, the convergence threshold is around p = 0.125. 

Fig.[8]shows two decoding trajectories superimposed on the 
EXIT chart of Fig. [7] at p = 0.125. We have used a 30, 000- 
qubit long interleaver. As seen from Fig. [8] the trajectory 


PTOI R (Outer) - PTOI R (Inner) 



Fig. 8. The EXIT chart of a QTC with decoding trajectories at p = 0.125 
(rate-1/9 QTC having PTOIR as both the inner and outer components with 
an interleaver length of 30, 000 qubits was used). 



Fig. 9. QBER performance curve with increasing iteration number for an 
interleaver length of 3,000 qubits. Rate-1/9 QTC having PTOIR as both the 
inner and outer components was used. 


successfully reaches the (x,y) = (1 ,y) point of the EXIT 
chart. This in turn guarantees an infinitesimally low QBER at 
p = 0.125 for an interleaver of infinite length. 

We have further verified the validity of our EXIT chart 
predictions using QBER simulations. Fig. [9] shows the QBER 
performance curve for an interleaver length of 3,000 qubits. 
The performance improves upon increasing the number of it¬ 
erations. More specifically, the turbo-cliff region starts around 
p = 0.125, whereby the QBER drops as the iterations proceed. 
Therefore, our EXIT chart predictions closely follow the 
Monte-Carlo simulation results. 

B. Entanglement-Assisted and Unassisted Inner Codes 

All non-catastrophic convolutional codes are non¬ 
recursive 0. Therefore, the resultant families of QTCs have 










































































a bounded minimum distance and do not have a true iterative 
threshold. To circumvent this limitation of QTCs, Wilde et 
al. ED, HD proposed to employ entanglement-assisted inner 
codes, which are recursive as well as non-catastrophic. The 
resulting families of entanglement-assisted QTCs have an 
unbounded minimum distance ina, ed, i.e. their minimum 
distance increases almost linearly with the interleaver length. 
Here, we verify this by analyzing the inner decoder’s 
EXIT curves for both the unassisted (non-recursive) and 
entanglement-assisted (recursive) inner convolutional codes. 

For classical recursive inner codes, the inner decoder’s EXIT 
curve reaches the {x,y) = (1,1) poin{§, which guarantees 
perfect decoding convergence to a vanishingly low QBER as 
well as having an unbounded minimum distance for the infinite 
family of QTCs (9) based on these inner codes. Consequently, 
the resulting families of QTCs have unbounded minimum dis¬ 
tance and hence an arbitrarily low QBER can be achieved for 
an infinitely long interleaver. This also holds true for recursive 
quantum convolutional codes, as shown in Fig. [10] In this fig¬ 
ure, we compare the inner decoder’s EXIT curves of both the 
unassisted and the entanglement-assisted QCCs of iflOl . which 
are labeled “PTOIR” and “PTOIREA”, respectively. For the 
PTOIR configuration, decreasing the depolarizing probability 
from p = 0.14 to p = 0.12 shifts the inner decoder’s EXIT 
curve upwards and towards the (1,1) point. Hence, the EXIT 
curve will manage to reach the (1,1) point only at very low 
values of depolarizing probability. By contrast, the EXIT curve 
of PTOIREA always terminates at (1,1), regardless of the 
value of p. Therefore, provided an open EXIT tunnel exists 
and the interleaver length is sufficiently long, the decoding 
trajectories of an entanglement-assisted QTC will always reach 
the (1,1) point; thus, guaranteeing an arbitrarily low QBER 
for the infinite family of QTCs based on these inner codes. 
In other words, the performance improves upon increasing the 
interleaver length; thus, implying that the minimum distance 
increases upon increasing the interleaver length and therefore 
the resultant QTCs have an unbounded minimum distance. 

C. Optimized Quantum Turbo Code Design 

The QTC design of fTol , lUTTl characterized in Fig. [8] 
exhibits a large area between the inner and outer decoder’s 
EXIT curves. The larger the ‘open-tunnel’ area, the farther 
the QBER performance curve from the achievable capacity 
limit (25). Consequently, various distance spectra based QTCs 
investigated in ED operate within 0.9 dB of the hashing 
bound. For the sake of achieving a near-capacity performance, 
we minimize the area between the inner and outer EXIT 
curves, so that a narrow, but still marginally open tunnel exists 
at the highest possible depolarizing probability. Our aim was to 
construct a rate-1/9 QTC relying on an entanglement-assisted 
inner code (recursive and non-catastrophic) and an unassisted 

5 Note that we only need (x. y) = (l,j/) for achieving decoding conver¬ 
gence to an infinitesimally low QBER. However, this requires an outer code 
having a sufficiently large minimum distance for the sake of ensuring that 
the outer code’s EXIT curve does not intersect with that of the inner code 
before reaching the (1, y) point. Unfortunately, an outer code having a large 
minimum distance would result in an EXIT curve having in a large open- 
tunnel area. Thus, it will operate far from the capacity. 



Fig. 10. Comparison of the inner EXIT curves of both unassisted and 
entanglement-assisted QCCs, labeled as PTOIR and PTOIREA respectively. 


outer code (non-catastrophic) having a memory of 3 and a rate 
of 1 /3. The resultant QTC has an entanglement consumption 
rate of 6/9, for which the corresponding maximum tolerable 
depolarizing probability was shown to be p max = 0.3779 

in ED- 

For the sake of designing a near-capacity QTC operating 
close to the capacity limit of 1X = 0.3779, we randomly 
selected both inner and outer encoders from the Clifford group 
according to the algorithm of (35) in order to find the inner 
and outer components, which minimize the area between the 
corresponding EXIT curves. Based on this design criterion, 
we found optimal inner and outer code pair whose seed 
transform^ (decimal representation) are given by: 

(Tinner = {4091,3736,2097,1336,1601,279, 

3093,502,1792,3020,226,1100}; (28) 

(Touter = {1048,3872,3485,2054,983,3164, 

3145,1824,987,3282,2505,1984}. (29) 

Fig. |TT] shows the corresponding EXIT chart at the conver¬ 
gence threshold of p = 0.35. As observed in Fig. QT| a 
marginally open EXIT tunnel exists between the two curves, 
which facilitates for the decoding trajectories to reach the 
(1,1) point. Hence, our optimized QTC has a convergence 
threshold of p = 0.35, which is only [10 x log 10 (o%ffg)] = 
0.3 dB from the maximum tolerable depolarizing probability 
of 0.3779. The corresponding QBER performance curves 
recorded for our optimized design are given in Fig. QT] A 
maximum of 15 iterations were used, while the interleaver 
length was increased from 1500 to 12,000. Similar to classical 
turbo codes, increasing the interleaver length for p < 0.35 
improves the attainable performance. Furthermore, Fig.fl2lalso 
compares our optimized design with the rate-1/9 QTC of ET1 

6 Please refer to on for the details of this representation. 
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Fig. 11. EXIT chart of the optimized rate-1/9 QTC (Interleaver length = 
30, 000 qubits). 



Fig. 12. QBER performance curves of the optimized rate-1/9 QTC for 
varying interleaver lengths and a maximum of 15 iterations. 


for an interleaver length of 3000, which is labeled “PTOIREA- 
PTOIR” in the figure. For the “PTOIREA-PTOIR” configu¬ 
ration, the turbo cliff region emerges around 0.31, which is 
within 0.9 dB of the capacity limit. Therefore, our EXIT-chart 
based QTC outperforms the QTC design based on the distance 
spectrum fTTl . More specifically, the “PTOIREA-PTOIR” 
configuration yields a QBER of 10 -3 at p = 0.2925, while 
our optimized QTC gives a QBER of 10 -3 at p = 0.3275. 
Hence, our optimized QTC outperforms the ‘PTOIREA- 
PTOIR” configuration by about [10 x log 10 ({y§§f§)] = 0.5 
dB at a QBER of 10~ 3 . However, our main design objective 
was to not to carry out an exhaustive code search, but to 
demonstrate the explicit benefit of our EXIT-chart based 
approach in the context of quantum codes. It must also be 
observed in Fig. [12] that a relatively high error floor exists 
for our optimized design, which is gradually reduced upon 


increasing the interleaver length. This is because the outer 
code has a low minimum distance of only 3. Its truncated 
distance spectrum is as follows: 

D{x) =2x z + 19x 4 + 10&r 5 + 530a: 6 + 2882a: 7 + 14179a; 8 + 
62288a; 9 + 243234a: 10 + 845863a' 11 + 1165784a: 12 + 
2501507a 13 + 744394a 14 . 

By contrast, the truncated distance spectrum of “PTOIR”, 
which has a minimum distance of 5, is given by DU: 

D{x) =lla 5 + 47a 6 + 253a 7 + 1187a® + 6024a 9 + 
30529a 10 + 153051a 11 + 771650a 12 . 

Consequently, as gleaned from Fig. [12] the “PTOIREA- 
PTOIR” configuration has a much lower error floor (< 10 6 ), 
since the outer code “PTOIR” has a higher minimum distance. 
However, this enlarges the area between the inner and outer 
decoder’s EXIT curves; thus, driving the performance farther 
away from the achievable capacity, as depicted in Fig. [8] 
Hence, there is a trade-off between the minimization of the 
error floor and achieving a near-capacity performance. More 
specifically, while the distance-spectrum based design primar¬ 
ily aims for achieving a lower error floor, the EXIT-chart based 
design strives for achieving a near-capacity performance. 

V. Conclusions 

In this contribution, we have extended the application of 
classical non-binary EXIT charts to the circuit-based syndrome 
decoder of quantum turbo codes, in order to facilitate the 
EXIT-chart based design of QTCs. We have verified the 
accuracy of our EXIT chart generation approach by comparing 
the convergence threshold predicted by the EXIT chart to the 
Monte-Carlo simulation results. Furthermore, we have shown 
with the aid of EXIT charts that entanglement-assisted recur¬ 
sive QCCs have an unbounded minimum distance. Moreover, 
we have designed an optimal entanglement-assisted QTC using 
EXIT charts, which outperforms the distance spectra based 
QTC of mi by about 0.5 dB at a QBER of 10~ 3 . 
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